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Abstract 

Using the a — uj model as an effective tool, the effects of self- 
consistency are studied in some detail. A coupled set of Dyson- 
Schwinger equations for the renormalized baryon and meson prop- 
agators in the a — uj model is solved self-consistently according to the 
dressed Hartree-Fock scheme, where the hadron propagators in both 
the baryon and meson self-energies are required to also satisfy this 
coupled set of equations. It is found that the self-consistency affects 
the baryon spectral function noticeably, if only the interaction with 
a mesons is considered. However, there is a cancellation between the 
effects due to the a and u mesons and the additional contribution 
of uj mesons makes the above effect insignificant. In both the a and 
a — uj cases the effects of self-consistency on meson spectral function 
are perceptible, but they can nevertheless be taken account of with- 
out a self-consistent calculation. Our study indicates that to include 
the meson propagators in the self-consistency requirement is unnec- 
essary and one can stop at an early step of an iteration procedure to 
obtain a good approximation to the fully self-consistent results of all 
the hadron propagators in the model, if an appropriate initial input 
is chosen. Vertex corrections and their effects on ghost poles are also 
studied. 
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1 . Introduction 



For a strong coupling theory it is important to go beyond the perturbation 
calculation. The Dyson-Schwinger (DS) equation provides one of the effective 
means to achieve such a goal. As is wellknown, the requirement of self- 
consistency can make the infinite series of Feynman diagrams contained in 
the solution to the equation nested, with each nest again containing the 
same infinite series of nested diagrams. However, the equation will then 
become nonlinear and not easy to solve. There is always the question what 
new effects on the solution may be produced by such a nesting and whether 
the latter is always important. In the zero-density case Brown, Puff and 
Wilets [1] as well as Bielajew and Serot [2] have solved the DS equation 
for the renormalized baryon propagators self-consistently according to the 
dressed Hartree-Fock (HF) scheme. As is wellknown, the a — uo model has 
achieved considerable success in the description of many bulk and single 
particle properties of nuclei [3]. Thus it is worthwhile to use it as a tool to 
study the effect of self-consistency in some detail. For simplicity we shall 
only consider the case of zero-density, where the tadpole self-energies need 
not be taken into account. The different self-consistent (dressed) HF schemes 
are illustrated in Fig. 1. By "dressed" it is meant that in the baryon and 
meson self-energies the hadron propagators are not all bare. If one requires 
self-consistency in the baryon propagators, i.e. G = G in Fig. la, while 
sets A = A and D = D , it is referred to as scheme BP. In this case the 
DS equation represented by Fig. la is already closed. Since G — G, it is 
easily seen by iteration that the self-energy is now nested and the nesting 
repeats again and again. Clearly one may also require self-consistency in 
all the hadron propagators contained in the self-energies. This implies that 
one should have G = G, A = A and D = D and it will be referred to as 
the fully self-consistent (FSC) scheme. In order to solve for G, A and D 
we must now consider a coupled set of DS equations, since the propagators 
are mutually correlated as shown in Fig. 1. Such a set of equations has 
been solved recently by Bracco, Eiras, Krein and Wilets (BEKW) [4] for a 
7i — uo model. Their calculation shows that the effect of self-consistency on 
the spectral properties of the hadron propagators is not significant and it also 
can not make the ghost poles disappear. Following previous works by Milana 
[5], Allendes and Serot [6] and Krein, Nielsen, Puff and Wilets [7], BEKW 
further found that the ghost poles in all the hadron propagators considered 
in the model can be eliminated if vertex corrections are properly taken into 
account. Recently, by means of a simple model we have shown [8] that the 
self-consistency in scheme BP diminishes the continuum part of the spectral 
representation for the baryon propagator and for the calculation of the latter, 
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there is a simpler scheme which is a good approximation to scheme BP. We 
would like to ascertain whether the above results hold more generally. We 
shall study the coupled set of DS equations according to the FSC scheme. 
The effects of self-consistency will be investigated in some detail. It is found 
that they affect the baryon spectral function noticeably, in case only the 
a meson is considered. However, there is almost no difference between the 
results obtained from scheme BP and FSC. It is interesting to find that 
the additional contribution of uj mesons makes the above effect insignificant. 
This indicates that there is a cancellation between the contributions from 
the a and uj mesons. In both the a and a — uj cases the effects of self- 
consistency on meson spectral functions are perceptible, but they can be 
taken account of without a self-consistent calculation. Our calculation shows 
that to include the meson propagators in the self-consistent requirement is 
unnecessary and one can stop at an early step of an iteration procedure to 
obtain a good approximation to the FSC scheme for the calculation of all 
the hadron propagators in the model. In analogy to Ref. [4], we have also 
obtained that neither scheme BP nor scheme FSC can eliminate the ghost 
poles. Since the question of ghost poles is of vital theoretical importance, 
following Refs. [4, 6, 7, 9], we have further calculated the contributions 
of vertex corrections and studied their consequences. As both the uj and 
a mesons are taken into account, vertex functions are only expressed by 
phenomenological form factors. 

In Section 2 we shall consider the coupled set of DS equations for the 
renormalized hadron propagators in the a — uj model and its reduction. The 
numerical results are given and discussed in Section 3. The vertex corrections 
and their consequences are studied in Section 4. Concluding remarks are 
presented in Section 5. 



2 . The model and coupled Dyson-Schwinger 
equations 

The Lagrangian density for the a — uj model [3] is given by 

L = -^{l,d, + M)V - \{d»4>d^ + mltf) - \F^ V 
-\m 2 v A^ + g s ipil)(j) + igyifrf^A^ + L CT c 

where F^ v = 0^A U U A^, 0^ = ~q^i = (•^j^o)) % = •^fj.^n = % %q 
with xq = t and CTC means the counterterm correction introduced for the 
purpose of renormalization. The baryon, a and uj meson propagators are 
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designated as follows: 

G a ?(x = Xl - x 2 ) = m a ( Xl )^(x 2 )]) = J ^e ik ^G a ?(k), (2) 

A(x) = (Tl^xMxi)]) = J ^e ik ^A(k), (3) 

D, u (x) = (T[A,( Xl )A u (x 2 )]} = J J^e^'D^k), (4) 

where (0) = (^ |0|^o) an d |^o) denotes the exact ground state. We note 
that our definition differs from the usual one by a factor of i. It is wellknown 
that the relevant DS equations in the dressed HF scheme (Fig. 1) can be 
written in the following form: 

(a) for baryon 

G{k) = G°(k) + G°(k)X(k)G(k) = -[7^ — %M + Z(k)]- 1 , (5) 
E(k)=E a (k) + E v {k), (6a) 
Z s (k) = -g 2 s J ^G(q)A s (k - q )F s (k, q,k-q) + ^ CTC (k), (6b) 

= 9l J -^l v G(q)D vX (k - q)T x (k, q,k-q) + ^ v CTC {k); (6c) 

(b) for a-meson 

A s {k) = A° s (k) + A^(k)U s (k)A s (k) = -i[k 2 + m 2 s + tU s (k) - ie\~\ (7) 

n s (k) = g] J ^~Jr[G(k + q)T s (k + q, q, k)G(q)} + U s CTC (k); (8) 

(c) for u;-meson 

= D%{k) + Dl v (k)U vX (k)D Xu (k), (9) 

%x(k) = -glj ^Tr[ lv G(k + q)T x (k + q,q,k)G(q)y, (10) 

where T s and T x are the a-baryon and o;-baryon vertex functions, respectively 
(denoted by a heavy dot in Fig. 1), a caret indicates that the quantity is 
not yet renormalized, r = 4 — 5 (5 — > + ) and in Eq. (5) the Feynman 
prescription M — > M — ie is understood. In this section vertex corrections 
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will not be considered, thus T s (p,q,k) = 1 and T\(p,q,k) = 7 A . Since the 
baryon current is conserved, we have 

n^(fc) = (^-^)n B (fc), (n) 

which implies II,, (fc) = \ H w (k). From Eq. (11) one observes that the 
renormalized H IJ , l/ (k) can be obtained from 

n v (k) = fi v (k) + ii v CTC (k). (12) 

To fix the renormalization counterterms we shall use the following conditions 



d 

s « ( k ) U^=iM t = 0; g( 7 k ) S "^^ 7 M fc M=^ = °> ( 13a ) 



n Q (k) \ k 2 =0 = 0; ^^IT a (A;)| fc2=0 = 0, (13b) 

where a = s or v and M t denotes the true baryon mass. If we set G = G, 
A s = A s and D v \ = D v \ in Eqs. (6, 8 and 10), it is seen that Eqs. (5- 
10) constitute a self-consistent closed set of DS equations for all the hadron 
propagators considered in the model. It thus yields the FSC scheme. We 
note that D vX can be written as [10, 11] 

D vX (k) = (S vX - -^) k2 + m i + m ^ k) _ ie - VkV^)' (14) 

where 5ml is the mass counterterm for the cu-meson. Since the u;-meson 
couples to the conserved baryon current, we may set D v \ = 5 V \A V in Eq. 
(6c), where according to Eq. (14) A„ is given by 

A v (k) = A° v (k) + A° v (k)U v (k)A v (k) = -i[k 2 +m 2 v + iU v (k) - ie]-\ (15) 

If we put A s = A° and D v \ = D° rjX and only require G = G, Eqs. (5) and 
(6) already build a closed self-consistent set, which has been called scheme 
BP. It is known that the original HF approximation to G(k) is given by 

G°M = -h^-iM + E(k,E k )}-\ (16) 

and the corresponding eigenvalue equation has the form 

[7^ -iM + Y,{k)] k0=Ek u{ks) = 0. (17) 
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If in addition to the approximations A s = A° and D v \ = 5 V \A®, we further 
substitute for G in Eq. (6), we obtain 

m = I (0 h ^ G ° s(g)A ° (A; - q)+ ^GUqh.Kik - ?)] + ^CTc(k), 

(18) 

Setting ko = E k in Eq. (18), one obtains the self-consistent equation for the 
HF potential Y*(k,Ek) = £(A;)|fc 0= £ fc [12]. Let the self-energy determined by 
Eq. (18) be indicated by a subscript p. We note that G p (k) = — — 

iM + E p (/c)] _1 is different from G^(k), because S P (A;) is defined for all k 
and k 0: though T, p (k)\ ko=Ek = £(£, E k ). It has been shown in Ref. [8] that 
G p (k) is a better approximation to G(k) determined by scheme BP than 
G^k) and such an approximation has been denoted as scheme P. Using G p 
we shall show that there is an extended potential scheme which yields a 
good approximation to the FSC scheme for the calculation of all the hadron 
propagators in the model. The four-dimensional integrals involved in Eqs. (6, 
8, 10) can be calculated by means of the spectral representation [1]. Following 
Ref. [2], we shall apply Feynman's integral parameterization to simplify 
the integration and use the standard method employed in the dimensional 
regularization to evaluate the relevant integrals. As is wellknown [13, 11], a 
properly normalized spectral representation for the baryon propagator may 
be expressed in the form 

r, n , v + 7 2 7 M V*(-m 2 ,Z 2 ) + iM t [3(-m 2 , Z 2 ) 

G(k) = -Z 2 j f 2 - / dm — ; , 

k 2 + Mr —ie J k 2 + m 2 — le 

(19a) 

where m± = M t + m s and Z 2 satisfies 

roo 

Z 2 + dm 2 a(-m 2 ,Z 2 ) = 1; (19b) 
whereas for the meson propagators we have 

Mk) = Z x , 2 ~l . - i I dm ^}- m ' Z *\ (20a) 
k 2 + mi — ie J k 2 + m 2 — ie 

th x 

where rh\ is the true meson mass, A = s or v, th\ means the threshold where 
the continuum starts and 

/•oo 

Z x + / dm 2 p x (-m 2 ,Z x ) = 1. (20b) 

Jth x 

Since in Eq. (1) we have only considered the linear meson interaction and no 
meson self-interactions are included, from Eqs. (8) and (10) it is easily seen 
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that th\ = (2M t ) 2 . Comparing Eq. (19a) with Eq._(5) and Eq. (20a) with 
Eq. (7) as well as (15), we have G(k) = ZG(k) and A x (k) = Y x A x (k), where 
Z and Y x are proportional factors. It is not difficult to see that Z 2 = ZZ t 
where (—Z t ) is the residue of G(k) at the pole 7„& iU = iM t , while under 
the on-shell renormalization condition (13a) one has Z t = 1 and M = M t 
[8]. Similarly, for the meson propagators Z\ = Y X R X , where (—iR\) is the 
residue of A x (k) at the pole k 2 = —rh\. Denote the denominator of A x {k) by 
D x (k 2 ) (U\ is a function of only k 2 ). The pole (— m 2 ) should clearly satisfy 
D x (-ml) = 0. From D x (k 2 ) = D x {-m\) + (k 2 + m 2 )^}^-* + • • •, one 
finds immediately 

R x = [i + i(dn x /dk?)„ = ^ l ]- 1 . (21) 

It is known [2] that in the zero-density case one may write 

S(jb) = -f^a(k 2 ) - iMb(k 2 ). (22) 
Denote a = a s + a v and b = b s + b v . Substituting 

G«c) = Z-G(t) = - 'lda ^ kM -f ) \ iM 'M-" 2 \ (23a) 



/ 7 (-a 2 ) = 6(a 2 - M 2 ) + 9(a 2 - m 2 ) 7 (-a 2 ), ( 7 = «or 0) (23b) 

A x (k) = Y x ^A x (k) = -ij dm 2 ^^- , (23c) 

o 

h x (-m 2 ) = R x 5{m 2 - m 2 ) + 0(m 2 - th x )p x (-m 2 ), (A = s or v) (23d) 

into Eq. (6) and following the procedure suggested in Ref. [2], one easily 
finds 

a s (A; 2 ) = tJ^// ^™ 2 / ^/ a (-a 2 )^(-m 2 )xln J + c s , (24a) 

oo ^ ' 

bs(k 2 ) = A- 2 Jj da 2 dm 2 J dxf,(-a 2 )h s (-m 2 ) In + C, (24b) 

oo \ t ) 

c - -Am 2 1 I da 2 dm 2 ) dr ^ - x) ^~ a ^ ± fp{-^)]h 8 {-m 2 ) 

° s ~ 8vr 2 ' JJ J K\-M 2 ) 

o o v ' 

(24c) 

K 2 (k 2 ) = K 2 {x, a 2 , m 2 , k 2 ) = x(l - x)k 2 + (1 - x)a 2 + xm 2 ; (24d) 
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a v (k 2 ) = 



b v {k 2 ) 



9v 
8tt 2 



4tt 2 



da 2 dm 2 / dxf a (—a 2 )h v (—m 2 )x\n 



K 2 (-M 2 ) 
K 2 (k 2 ) 



+ c„, (25a) 



CXJ 1 

J J da 2 dm 2 J dxfp{—a 2 )h v {—m 2 )hi 



K 2 (-M 2 ) , x 



4tt 2 



dcr <im / dx 



x(l - x)[2fp(-a 2 ) - x f a (-a 2 )}h v (-m 2 ) 



K 2 {-M 2 ) 



where ^-a 2 ) = Z'^-a 2 , Z 2 ), Px (-m 2 ) = Y^ 1 p x (-m 2 , Z \) and 



a(k 2 ) = —Im 

7T 



(3{k 2 ) = -Im 



l + a s (k 2 ) + a v (k 2 ) 
D(k 2 ) 

l + b s (k 2 ) + b v (k 2 ) 
D(k 2 ) 



(25c) 

(26a) 

(26b) 
(26c) 



D{k 2 ) = (1 + a s (k 2 ) + a v (k 2 )) 2 k 2 + (1 + b s (k 2 ) + b v (k 2 )) 2 M 2 . 

It is seen that h\(—m 2 ) is yet unknown and use must now be made of Eqs. 
(7-10). Since Eqs. (8) and (10) can be worked out in the same way, let us 
consider the latter. Substituting Eqs. (23a and b) into Eq. (10) and using 
the Feynman formula (AB)~ l = f£° dx[xA + (1 — x)B]~ 2 , one gets 



E„n™(fc) = 8^ 



da 2 dm 2 I dx I 



\{Q 2 - x(l - x)k 2 )F a {a 2 m 2 ) + 2M 2 F (3 (a 2 m 2 )}[Q 2 + K 2 {k 2 )\ 



(27) 



where F 7 (a 2 m 2 ) = / 7 (— <r 2 )/ 7 (— m 2 ). By means of Eq. (13b) and the inte- 
gration formula familiar in the dimensional regularization we obtain that 



n,(& 2 ) = lE v iK(k 2 ) - n w (o) - k 2 (dfi v jdk 2 ) k2=0 ] 

= i£s IT ' da 2 dm 2 f dx{[(K 2 (k 2 ) + \x{\ - x)k 2 )F a (a 2 m 2 )- (28) 
M 2 F,(a 2 m 2 )} lnf@ - x(l - x)k 2 [F a (a 2 m 2 ) - ^F,(a 2 m 2 )]}, 

where the counterterms are supplied by terms |5m 2 A M A M and <, V \F^ U F^ V 
contained in L C tc [ see Eq. (1)]. Similarly we have 

T[ s (k 2 )=i^ jr da 2 dm 2 ^ dx{[(2K 2 (k 2 ) + x{\ - x)k 2 )F a (a 2 m 2 ) + 
M 2 F^a 2 m 2 )} lnf^ - x(l - x)k 2 [2F a (a 2 m 2 ) + jgfeF^m 2 )]}. 

(29) 
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Further, using Eqs. (7, 15, 20) and the formula (x — ie) 1 = f + irr5(x), we 
obtain 

p a (k 2 ) = -Im[iA a {k% (30) 

7T 

for k 2 < —th\. Eqs. (24-26) and (28-30) are the explicit expressions for the 
closed set of renormalized DS equations used for the calculation. We note 
that through the renormalization conditions (13a and b) G(k) and A a (k) 
given by Eqs. (5, 7, 15) are normalized differently from G(k) and A\(k), 
respectively. According to Eqs. (19b) and (20b) we shall define h (I x ) given 
below _ 




h = Y x dm 2 Px {-m 2 ) = 1 - Y X R X , (20c) 

Jth x 

as the impurity factor of the single particle (sp) character of a baryon (A- 
meson). Clearly, it just tells the relative importance of the continuum part. 

3 . Numerical results 

Though the set of self-consistent DS equations for scheme FSC obtained in 
Section 2 looks somewhat formidable, it can be solved rigorously and quite 
easily by the method of iteration. Thus, it provides a useful example for 
the study of self-consistency and diagram nesting. In the following the mass 
values used are M t = 4.7585, m s = 2.6353 and m v = 3.9680(/m- 1 ). We 
shall first consider the case of baryons interacting only with a mesons. Two 

2 

values of the strength parameter g 2 (= yf^) are studied. Since the results 
of ~g 2 = 0.5263 and 0.6517 display the same qualitative behavior, only the 
former will be represented graphically. From Figs. 2 and 3 one sees that 
the results obtained from schemes BP and FSC are almost the same. Let us 
denote a physical quantity Q(k) by Q(k; S) or simply Q(S), if it is calculated 
according to scheme S. Though one has set A s = A° in scheme BP, one may 
substitute the baryon propagator G(k; BP) obtained from it into Eq. (8) and 
find p s (k 2 ; BP) from Eq. (30). It is found that p s (k 2 ; BP) and p s (k 2 ; FSC) 
are very close to each other. Thus, only the latter is shown in Fig. 4. The 
above results indicate that to include the meson propagators in the self- 
consistency requirement may not be necessary. Form Eqs. (5-10) it is seen 
that all the self-energies will be known, if G, A s and D vX are known. Thus 
different approximate schemes are obtained, if different choices are made for 
G, A s and D r/X . To assess the importance of self-consistency and to ask 
whether one may stop at an earlier step of the iteration process, we have 
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shown in Fig. 3 results from five different calculation schemes which are 
explained in Table 1. The first column gives the name of each scheme, while 
the second and third explain how its S s [E„] and IT S [ITJ are obtained from 
Eqs. (6b) [(6c)] and (8) [(10)], respectively. For instance, S S (P) and H S (P) 
in scheme P are obtained by setting G = G^, A s = A° in Eq. (6b) and 
G = Gl in Eq. (8), Z s (PEP)_and Tl s (PEP) in the partially extended P 
scheme (PEP) by G = G(P), A s = A° s in Eq. (6b) and G = G(P) in Eq. 
(8), whereas to obtain Tj s (EP) and U S (EP) for the extended P scheme (EP), 
one sets G = G(P), A s = A S (P) in Eq. (6b) and G = G(P) in Eq. (8), 
etc. The quantities in the brackets are for cu-mesons in the a — uo model. 
From Fig. 3 it is seen that the self-consistency diminishes the continuum 
part of the spectral representation for the baryon propagator [8] , though the 
convergent process from (a(P), (3{P)) to (a(FSC), /3(FSC)) is oscillatory. 
It is interesting to note that the contribution of the meson self-consistency, 
though very small, adds toward the same direction. Clearly we have 

a(PEP) ~ a(EP), (3{PEP) ~ (3(EP) 
a(BP)~a(FSC), P(BP) ~ P(FSC) ' { } 

whereas Fig. 4 shows 

Ps (PEP) = Ps (EP) ~ Ps {BP) ~ Ps (FSC). (31b) 

In order to gain some quantitative insight into the above relations we have 
calculated and listed the values (J CT ) of the baryon (a-meson) sp impurity 
factor in Tables 2 and 3. In fact, the small values of I a also explain why it 
is not important to require self-consistency in the meson propagators. One 
observes that the difference between the results of schemes PEP and EP is 
small and PEP seems a quite good approximation to FSC with respect to 
all the hadron propagators considered. To ascertain this we have further 
calculated 

F(±) = (iG(k))± = ut{ks)iG{k)u±{ks) 



l+q(fc 2 ) 



l+a(k 2 ) - 1 - E k 



M t < 32 > 



D(fc 2 ) 

[8], where +(— ) refers to the eigenspinor of Eq. (17) with eigenvalue +(— )E^ = 
[M?+k 2 ]^. In Fig. 5 we have only plotted F r ^(FSC) = Re(Im)(iG(FSC)} + , 
with k = v / 350/m~ 1 , since the numerical results of the three schemes P, 
PEP and FSC are nearly the same. Thus, if we use Eq. (32) as a criterion, 
we note that even G(P) is a good approximation. The poles of G(k) and 
A\(k) are given by the roots k 2 of D(k 2 ) of Eq. (26c) with a v = b v = and 

F > x(k 2 ) = k 2 + m\ + ill\(k 2 ), (33) 
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respectively. The ghost poles found are listed in Table 4. They are given in 
units of fmT 2 . 

Now let us consider the a — uj model. Here we shall only present the 

2 

numerical results for g 2 = 0.5263 and g 2 , = = 1.3685, as qualitatively 
they are quite typical. In Figs. 6 and 7 we have plotted the numerical results 
for (a s ,b s ) and (a v ,b v ). It is seen that before k 2 ~ — 120/m~ 2 the effects of 
self-consistency is insignificant. Hereafter they seem to become more impor- 
tant, but they are not large enough to cause any palpable effects on a and (5. 
The results obtained from the five schemes nearly coincide with each other 
as shown in Fig. 8, where only two are drawn explicitly, because the other 
four are not distinguishable in the figure. Thus, in this case no self-consistent 
requirement is needed, just as found in Ref. [4] for the n — uj model. Compar- 
ing Figs. 6 and 7 with Fig. 2 one notes that a s (k 2 ) and b s (k 2 ) only change 
slightly due to the presence of uj mesons. However, there is a quite large 
difference between (a s , b s ) and (a, b) where a = a s + a v and b = b s + 6„.This 
explains why Fig. 8 differs from Fig. 3 so greatly. A comparison between 
Figs. 3 and 8 suggests that there exists a cancellation between the effects on 
self-consistency due to the a and uj mesons. This is obviously related to the 
a-baryon and u;-baryon interactions and thus also to the different character- 
istic features contributed to the baryon-baryon interaction by the a and uj 
mesons individually. It seems that these features are not only of importance 
to the nuclear binding and saturation, but also to the formation of baryonic 
excited states (see below) . We expect the above relation between the cancel- 
lation and the different characteristics of meson-baryon interactions may hold 
generally, though the details of its regularity will depend on the participating 
meson fields and there may even be enhancement if their characteristics are 
of the same kind (which is actually defined by the 'enhancement'). If the 
regularity of cancellation can be extended to the case of finite density, it will 
be very useful, because the analysis and the renormalized calculation will be 
greatly simplified, when the effect of self-consistency is negligible. However, 
this is still a problem to be studied. To illustrate how the uj meson affects 
the disappearance of the effect of self-consistency we have drawn a set of a 
and P calculated for an intermediate g 2 = 0.3400 in Fig. 9. We note that 
owing to the presence of uj mesons the maximum of each baryon spectral 
function becomes sharper and more distinct. Since a{k 2 ) relates directly to 
the probability of occurrence of an excited baryon state, it indicates that 
the contribution of uj mesons enhances the possibility of forming a resonant 
baryon state, though we do not expect that the a — uj model may explain 
the experimentally observed baryon spectrum. From Fig. 10 it looks that a 
discernible effect of self-consistency exists for mesons. However, we observe 



11 



from the same figure that schemes PEP and FSC yield almost the same result. 
This says that the above effect can be taken account of by simply substituting 
G(k; P) in Eqs. (8) and (10) without considering a self-consistent calcula- 
tion. BEKW [4] found that for the 7r — uj model the self-consistency decreases 
p n (k 2 ), while increases p u (k 2 ). Owing to the difference in models, our cal- 
culation shows that the self-consistency decreases both p s (k 2 ) and p v (k 2 ). 
Further we note that the self-consistency and additional contribution of uj 
mesons still cannot eliminate the ghost poles, whose values found for the rel- 
evant hadron propagators are respectively, for baryon: 149.4215 ± i96.8997, 
for a meson: 274.4833 and for uj meson: 178.5431 (fm~ 2 ). 



4 . Effects of vertex corrections on ghost 
poles 

The principal aim of this section is to study the effects of vertex corrections 
on ghost poles and their elimination. For our purpose we shall only consider 
phenomenological form factors, though they still suffer the not yet resolved 
problem of violating the Ward-Takahashi identity as emphasized in Refs. [3, 
7]. The vertices included in Eq. (1) are illustrated in Fig. 11, where 1 and 
2 denote baryons, while A a a or uj meson. Phenomenologically the vertex 
functions may be written approximately as [14] 

r s (pi,p 2 ,P s ) = F b (p 1 )F b (p 2 )F s (p s ), (34a) 

r„(pi,P2,P„) = l ll Fb{pi)F b {p2)F v (p v ), (34b) 

where F b in T s and T M will be chosen the same. Substituting Eq. (34) in 
Eqs. (6, 8, 10), we obtain 

Eftfc) = -9 2 s F b {k) J ^L Fb (q)F s (k - q)G(q)A s (k - q), (35a) 

Ef(fc) = g 2 v F b (k) J -^L Fb (q)F v (k - q) lfl G(q)D, u (k - q) lv , (35b) 

Hf(fc) = 9 2 s F s (k) J ^ 4 F b ( k + q)F b {q)Tr[G{k + q)G(q)}, (36) 

n^(fc) = ~9 2 v F v (k) J j^y 4 F b ( k + q)F b (q)Tr[ ltl G(k + q) lv G{q)\, (37) 
where the superscript / indicates that form factors are inserted. 
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In the previous section it has been shown that all the calculated renor- 
malized hadron propagators possess ghost poles. Just as found in Refs. [4, 
6, 7, 9], our calculation also shows that the ghost poles can be eliminated by 
a proper treatment of vertex corrections. In order to make the change of the 
/c-dependence of the self-energies caused by them appear in a clearer form, 
we shall choose the form factors as simple as possible. We shall set F b — 1 
and 



m 2 

= A2 i 2 A . , \ = s or v (38) 

where e — > + indicates the displacement of the poles. From Eq. (35) it 
is immediately seen that S{(A = s or v) are made finite by this choice. 
However, according to Eqs. (36) and (37) we have 

lift*) = F s (k)g 2 J ^Tr[G(k + q)G(q)] = F a (k)U a (k), (39) 

n>,(fc) = -F v {k)g 2 v J ^Tr[ 7/1 G(fc + q) lv G{q)\ = F v (k)U^(k) : (40) 

where U s (fl^) is just the integral in Eq. (8) (Eq. (10)) with T s = 1(T U = 
7„). Both of them are divergent. As no confusion will arise, the bar over 
will hereafter be omitted. In order to investigate the effect of the multiplica- 
tive factor F\ on the elimination of ghost poles, we shall, just as Ref. [6], 
substitute the renormalized IT S and IT M1 , for II S and fl^, respectively. The 
explicit expressions of IT^ and IJ S have been given in Eqs. (11, 12, 28, 29). 

Now let us consider the baryon self-energy. Substituting Eq. (38) into 
Eq. (35) and using the method of integral parameterization, we get 

Ef(k) = -2ig 2 JJda 2 dm 2 (A 2 -m 2 )h s (-m 2 ) 

d^Q (x+y^^fcj-a^+iMtM-a 2 ) , (41a) 



1 fl-X 



o J o V J ^ W+^j-in] 



L 2 = (x + y) (1 - x - y) k 2 + xA 2 s + ym 2 + (1 - x - y)a 2 , (41b) 

where 77 — > + . It is seen that the Wick rotation can be applied and the 
four-dimensional integral can be worked out immediately. From Eqs. (22) 
and (41) we obtain 



af s( k ) ^ jjl ,h:1, ' in ~ j J' (A^T 



xf a (-a 2 )h s (-m 2 ) In 



L 2 (z,z,a 2 ,m 2 ,k 2 ) ' ( 42& ) 



H(k) = if //J doHm 2 f Q dz 

xf p (-a 2 )h s (-m 2 ) In L l [z ^^^ 



(A|— m|) 

(42b) 
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L){x, z, a 2 , m 2 , k 2 ) = z(l - z)k 2 + (1 - z)a 2 + zm 2 + x(A 2 s - m 2 ). (42c) 

According to Eq. (42c) L 2 (z, z,a 2 ,m 2 ,k 2 ) is independent of m 2 and in the 
neigbourhood of m 2 = A 2 

L}(z, z, a 2 , m 2 , F) L 2 (^, z, a 2 , m 2 , F) l ™ J + " ' ' 1 } 

thus the integrand is well behaved at m 2 = A 2 even if A 2 > th s . Similarly 
from Eq. (41) 



al(k) = & fftdo*dm* \\dz ^ z 

w t I 2\L / 2 \ 1 L 2 ,(z,z,cr 2 ,m 2 ,k 2 ) ' 

x U-a )K{-m)\n ^ z(j2m ^ k2) 

H(k) = Hi JJ~do*dm* J\dz 

xffs(-a 2 )h v (-m 2 ) In jj^J^J^ 



(44a) 



(44b) 



Eqs. (42) and (44) combined with Eqs. (23, 26, and 28-30) build a closed 
set, which has been solved by iteration. However, one notes that instead of 
Eqs. (26b and c) one should now write 



l + b s (k 2 ) + b v (k 2 ) 



(26d) 



D(k 2 ) 

D(k 2 ) = (1 + a s (k 2 ) + a v (k 2 )) 2 k 2 + (1 + b s (k 2 ) + b v (k 2 )) 2 M 2 , (26e) 
while between M t and M we have 

(1 + a s (-M 2 ) + a v (-M 2 )) 2 M 2 = (1 + 6 S (-M 2 ) + b v (-M 2 )) 2 M 2 . (27) 

The numerical results are represented graphically in Fig. 12 to 15 for A s = 
A v = 5.0676/m~ 1 (lGel / ). The other parameters are the same as chosen 
previously. According to our calculation no noticeable difference between 
scheme BP and FSC exists, thus only the latter is plotted. From the fig- 
ures it is seen that the self-consistency is almost of no effect. Comparing 
Figs. 12-13 with Figs. 6-7 one notes that the two set of (a, b) differ signif- 
icantly, especially with respect to their asymptotic behavior. However, the 
results given in Figs. 8 (10) and 14 (15) for the baryon (meson) spectral 
functions are qualitatively similar. In Table 5 we have presented the values 
of lb and I\ calculated according to the method of renormalization and of 
form factors. If we adjust the parameters, it is possible to make the two 
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sets approach each other better. This suggests that for the calculation of 
hadron propagators the two methods can yield closely approximated results. 
We have found that all the ghost poles are eliminated. The reason seems 
intimately related with the change of the asymptotic behavior caused by 
the vertex correction. For the baryon self-energy we see that as \k 2 \ — > oo, 
the Log-function | \n[K 2 (- M 2 ) / K 2 (k 2 )}\ -> oo in Eqs. (24, 25), while in 
Eqs. (42, 44) we have | \n[L}(z, z, a 2 , m 2 , k 2 )/L}(0, z, a 2 , m 2 , k 2 )}\ -> 0. This 
fact also explains the asymptotic behavior mentioned above. As shown in 
Eqs. (39, 40) the meson self-energies Il{ are obtained by simply multiplying 
the corresponding expressions given in Section 2 by F\(k 2 ), which tends to 
(A; 2 !" 1 as \k 2 \ -> oo. From Fig. 16 and Eqs. (28, 29) one sees that iU x (k 2 ) 
bends the straight line k 2 + m\ to intersect the abscissa axis (/c 2 -axis) and 
thus A\(k) gets a ghost pole. However no such bending will be induced by 
iU. x = F\{k 2 )iTi\{k 2 ) as shown in Fig. 17, because if k 2 becomes large, 
will be negligible compared with k 2 . It is interesting to note that there will 
be no ghost poles if iU\(k 2 ) bends the straight line upward or the curve 
D\ = D\(k 2 ) turns upward before it intersects the /c 2 -axis. 

5 . Concluding remarks 

Under the self-consistent HF approximation we have solved the coupled set 
of DS equations for the renormalized hadron propagators in the a — oo model. 
Since this set can be solved rigorously, it provides a convenient means for the 
study of the effects of self-consistency. Our calculation shows that in the case 
of zero-density scheme PEP gives a quite good approximation to the FSC 
scheme and there is no need to require self-consistency in the meson propaga- 
tors. In the cr — uj case the self-consistency almost has no effect on the baryon 
propagator, since there is a cancellation between the effects caused by the 
two kinds of mesons, whereas for the a case a FSC calculation is also unnec- 
essary, because the iteration procedure converges quickly if an appropriate 
initial input is chosen. Fig. 3 shows that the self-consistency diminishes the 
continuum part of the baryon spectral representation predicted by scheme P. 
From Tab. 2 one observes that this part is by no means always unimportant. 
However, there is generally no need to consider scheme BP, since if scheme 
P does not yield a good enough result, one may go one step further to use 
scheme PEP. For the study of nuclear matter and nuclear structure one has 
to take the case of non-zero density into account. The question whether in 
the latter case the above conclusions still hold remains to be studied. Us- 
ing the method of phenomenological form factors, we have also investigated 
the significance of vertex corrections. Our calculation confirms the results of 
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Refs. [4, 6, 7, 9], namely the ghost poles in the hadron propagators can be 
eliminated by vertex functions with damping asymptotic behavior. In Ref. 
[7] it was found that the baryon spectral function Ar(k) can be negative for 
some values of real k if only the coupling with cj-mesons is considered. In the 
latter case we have also found that a(k 2 ) may be negative in a region close 
to the threshold. Though, if in addition to u other mesons are taken into 
account, the above problem does not occur, it still has to be studied and is 
being under study, because the coupling with u;-mesons is of importance and 
according to their definitions and physical meanings, both Ar(k) and a(k 2 ) 
should be > 0. 

The work is supported in part by the National Natural Science Foundation 
of China and the Foundation of Chinese Education Ministry. 
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Table 



Table 1 : Different calculation schemes. 



Name 




n s (n„) 


P 


G u s , (A°) 




PEP 


G(P), A° s (AS) 


G(P) 


EP 


G(P), A°(P) (A°(P)) 


G(P) 


BP 


G(PP), A° (A°) 


G(PP) 


FSC 


G(FSC), A°(FSC) (A°(FSC)) 


G(FSC) 



Table 2 : Values of the baryon impurity factor 



Scheme 


P 


PEP 


EP 


BP 


FSC 


^ = 0.5263 


0.4157 


0.3528 


0.3504 


0.3672 


0.3664 


^ = 0.6517 


0.4638 


0.3253 


0.3227 


0.3756 


0.3752 



Table 3 : Values of the a-meson impurity factor I a . 



Scheme 


P 


PEP 


FSC 


g'i = 0.5263 


0.1725 


0.1178 


0.1241 


ts = 0.6517 


0.1557 


0.0931 


0.1017 



Table 4 : Ghost poles of the baryon and a-meson propagators. 



Propagator 


Baryon 


cr-meson 


g\ = 0.5263 


35.5356 ± i253.5311 


246.7354 


g 2 s = 0.6517 


-13.6188 ±il27. 7365 


177.4338 



Table 5 : Impurity factors calculated according to the method of renormal- 
ization (I) and of form factors (II). 



Impurity factor 


h 


la 


L 


I 


0.2905 


0.1334 


0.1870 


II 


0.2701 


0.0459 


0.0643 



Table 1 on P 9 , Tables 2-4 on P 10 , Table 5 on P 14 . 
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Figure captions 

Fig. 1 : Diagrammatic representation of the different self-consistent (dressed) 
HF schemes, a. the baryon propagator, b. the cr-meson propagator, c. 
the a;-meson propagator. 

Fig. 2 : Numerical results of (a s , b s ). (a) the real part, (b) the imaginary 
part. [BP ~ FSC, ~ means indistinguishable or almost indistinguish- 
able.] 

Fig. 3 : The baryon spectral functions a(k 2 ) and (3(k 2 ). [a(BP) ~ a(FSC); 
a(PEP) ~ a(EP)\ 

Fig. 4 : The cr-meson spectral functions p s {k 2 ). 

Fig. 5 : Graphical representation for F(+) =< iG(k) > + , where F r = Re < 
iG HF (FSC) >+, F r = Im< iG HF (FSC) >+. 

Fig. 6 : Numerical results of a s and ) the real part, (b) the imagi- 

nary part. [a sr (BP) ~ a sr (FSC); a si (BP) ~ a si (FSC); a vr (BP) ~ 
a OT (F5C)] 

Fig. 7 : Numerical results of b s and 6„. (a) the real part, (b) the imag- 
inary part. [b sr (BP) ~ 6 sr (F,SC); & si (SP) ~ b si (FSC); b vr {BP) ~ 

Fig. 8 : The baryon spectral functions a(/c 2 ) and /9(/c 2 ) calculated with 
baryon coupling with both a and oj meson. 

Fig. 9 : The baryon spectral functions a{k 2 ) and j3(k 2 ) for g 2 = 0.34. 

Fig. 10 : The meson spectral functions p x (k 2 ). (a) cr-meson: \ = a; (b) 
cu-meson: A = u. [PEP ~ FSC] 

Fig. 11 : The baryon-meson vertices. 

Fig. 12 : Numerical results of (a s , a v ) [FF]. where the square bracket [FF] 
indicates that the results are obtained by the method of form factors. 

Fig. 13 : Numerical results of (b s , b v ) [FF]. 

Fig. 14 : The baryon spectral functions a(k 2 ) and (3(k 2 ). 

Fig. 15 : The meson spectral functions p x (k 2 ) [FF]. (a) cr-meson: A = s; 
(b) oj- meson: A = oj. 
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Fig. 16 : Plot of D\ = D\(k 2 ) versus k 2 . (a) cr-meson: A = s; (b) u;-meson: 
A = cu. 

Fig. 17 : The same as Fig. 16 but calculated by the method of form factors. 

Fig. 1 on P 2 , Figs. 2-4 on P 9 , Figs. 5-7 on P 10 , Figs. 8-10 on P n , Figs. 
12-17 on P M . 



20 



0.002 - 



0.001 - 



0.000 - 



0.001 - 



■0.002 




P(k 2 ) 



■600 



■500 



■400 



■300 



Fig. 




P-Scheme 
PEP-Scheme 
EP-Scheme 
BP-Scheme 
FSC 



■200 



100 | 2 fC -2v 

k (fin ) 



0.0004 

P c (k 2 ) 



0.0003 - 



0.0002 - 



0.0001 - 



0.0000 



1 1 ' 

:P-Scheme 


i • i • i • 

R =0.9374 


■ :PEP-Scheme 


R =0.9150 '\ 

s 


:FSC 


R =0.9183 

s 

• • 

• 






• • 

" s * \ 

/' V 

/' V. 
/' X 

/> V 

s } 

y ' 

— ^ ' 

. — ' 




i 


i.i.i 





1000 



■800 



■600 -400 

Fig.4 



- 200 k 2 (fm" 2 ) 



1.5 



1.0 - 



0.5 - 



0.0 



■0.5 



1 1 1 1 1 

■ 


i • i • i • i 1 

F 

r \ 




F 

. i 


F \ 

r \ 

i . i . i . i . i . i \ 



Fig. 5 



-350 -300 -250 -200 -150 -100 



■50, 2 /£ -2 X 

k (fm ) 



0.0030 - 



0.0025 - 



0.0020 - 



0.0015 - 



0.0010 - 



0.0005 - 



: P-Scheme 
: BP-Scheme 



0.0000 




■600 




:P-Scheme 
:PEP-Scheme 
:FSC 



R =0.9374 

s 

R =0.9212 

s 

R =0.9227 



0.0006 

P v (k 2 ) 

0.0004 



:P-Scheme 
:PEP-Scheme 
:FSC 



R =0.7772 

V 

R v =0.7424 
R =0.7449 



0.0002 



0.0000 



1000 



-800 



■600 -400 

Fig. 10 



(a) 





(b) 



0.0004 

• P s (k 2 ) 

- 0.0003 



- 0.0002 



- 0.0001 



0.0000 



-200 | 2 /r -2. 

k (fm ) 




Fig. 1 1 



0.006 



0.005 - 



0.004 - 



0.003 - 



0.002 - 



0.001 - 



0.000 




-400 



-300 



-200 

Fig. 14 



100 ,2 /r -2v 

k (fm ) 



0.00060 

P,(k 2 ) 



0.00045 - 



0.00030 - 



0.00015 - 



0.00000 



:P-Scheme 
:PEP-Scheme 
:FSC 



R=0.9624 
R=0.9525 

s 

R =0.9552 




:P-Scheme 
:PEP-Scheme 
:FSC 



R v =0.8004 
R =0.7696 

V 

R =0.7802 



i.iiii.. T* 8 * 




■600 



-500 



-400 



-300 

Fig. 15 



-200 



(a) 



0.00020 

Ps* 2 ) 

-10.00015 



- 0.00010 



- 0.00005 



(b) 



0.00000 



100 k 2 (fm 2 ) 



60 - 



-20 L 



' 1 1 1 1 1 1 1 1 1 1 

^_ (a) 
/ D =k 2 +m 2 +iIl (k 2 ) \ 

_ / s s s v 7 \ 


1 1 1 1 1 1 1 1 1- 

(b). 

D =k 2 +m 2 +in (k 2 ) \ 

V V V \ 


^\ ■ 

(178.5431,0) \ 



A 

(274.4833,0) ^ \ 




- 40 



- 20 



- -20 



- -40 



^-60 



50 100 150 200 250 300 50 100 150 200 

k 2 (fm 2 ) 



Fig. 16 




k 2 *10 4 (fm" 2 ) 
Fig. 17 



